Evolution Effects on Parton Energy Loss with Detailed Balance 
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The initial conditions in the chemical non-equilibrated medium and Bjorken expanding medium 
at RHIC are determined. With a set of rate equations describing the chemical equilibration of 
quarks and gluons based on perturbative QCD, we investigate the consequence for parton evolution 
at RHIC. With considering parton evolution, it is shown that the Debye screening mass and the 
inverse mean free-path of gluons reduce with increasing proper time in the QGP medium. The 
parton evolution affects the parton energy loss with detailed balance, both parton energy loss from 
stimulated emission in the chemical non-equilibrated expanding medium and in Bjorken expanding 
medium are linear dependent on the propagating distance rather than square dependent in the 
static medium. The energy absorption can not be neglected at intermediate jet energies and small 
propagating distance of the energetic parton in contrast with that it is important only at intermediate 
jet energy in the static medium. This will increase the energy and propagating distance dependence 
of the parton energy loss and will affect the shape of suppression of moderately high Pt hadron 
spectra. 

PACS numbers: 12.38.Mh,24.85.+p,25.75.-q 
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I. INTRODUCTION 



One of the challenging goals of heavy-ion physics is to 
detect quark-gluon plasma (QGP). The two nuclei pass 
through each other, interact, and then produce a dense 
plasma of quarks and gluons. As the initial parton den- 
sity is large and the partons suffer many collisions in a 
very short time, the initial partonic system may attain 
kinetic equilibrium. But does it attain chemical equilib- 
rium? This question has been investigated in the frame- 
work of parton cascade model Q ,which is based on the 
concept of inside-outside cascade [2|-|4[ and evolve parton 
distributions by Monte-Carlo simulation of a relativis- 
tic transport equation involving lowest order perturba- 
tive QCD scattering and parton fragmentations. From 
the numerical studies (5l4lCj| three distinct phases of par- 
ton evolution can be distinguished: (1) Gluon thermalize 
very rapidly, reach approximately isotropic momentum 
space distribution after a time of the order of 0.3 fm/c. 
(2) Full equilibration of gluon phase space density takes 
considerably longer. (3) The evolution of quark distribu- 
tions lags behind that of the gluons because the relevant 
QCD cross sections is suppressed by a factor of 2-3. This 
calculation indicates that the QGP likely to be formed 
in such collisions are far from chemical equilibrium. 

Gluon radiation induced by multiple scattering for an 
energetic parton propagating in a dense medium leads to 
induced parton energy loss or jet quenching. Jet quench- 
ing is manifested in both the suppression of single in- 
clusive hadron spectra at high transverse momentum px 
region and the disappearance of the typical back-to- 
back jet structure in dihadron correlation as discovered 
in high-energy heavy-ion collisions 
theoretical investigation of jet quenching has been widely 
carried out in recent years Qj§-[l8| . It is found that in the 



static medium the radiative energy loss is proportional to 
square of propagating distance. Later the detailed bal- 
ance effect with gluon absorption was included. It has 
been shown that the gluon absorption play an impor- 
tant role for intermediate jet energy region in the static 
medium fl9j . But the plasma is not static, it will expand, 
cool and become more dilute. Recently medium expan- 
sion is included in the jet energy loss by stimulated emis- 
sion in Ref. (20| . However, only the gluon distribution 
evolution in a thcrmodynamical equilibrated expanding 
medium is considered, the temperature evolution, which 
play an important role for debye screened mass calcu- 
lation and affects the mean free path and opacity, was 
neglected in Ref. [2(J. Moreover, since the QGP is like 
to be far from chemical equilibrium, the effect of par- 
ton chemical equilibration on jet energy loss need to be 
studied. 

In this paper, we will study the effects of tempera- 
ture and fugacity (chemical potential) evolution on jet 
quenching with detailed balance. In general, the question 
of thermodynamical equilibration can be decided with 
microscopic transport models [5]. In this paper, we study 
a simpler problem which allows us to use a macroscopic 
model to get the initial conditions and parton evolution. 
Our strategy will be to assume that the parton distribu- 
tion can be approximately by thermal phase space distri- 
bution with non-equilibrium fugacity X g and X q of gluon 
and quark, use a set of rate equations to describe the 
chemical equilibration of partons, compare dEx/dy and 
dN/dy which we get with the data from RHIC, we can 
then determine the initial conditions of proper time, fu- 
gacity, and the time dependence of the parameters T, \ g 
at RHIC [12|. The and X q . We will also study the a thermal-equilibrated 



and chemical-equilibrated medium expansion - Bjorken 
expansion with \ g ( q ) = 1 and compare the difference. 
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With these results, we obtain the Debye screening mass, 
mean free path, and opacity from the perturbative QCD 
at finite temperature in a thermal equilibrated, but chem- 
ical non-equilibrated medium and find that they are dif- 
ferent with those in thermodynamical medium and static 
medium. Then, we will investigate the evolution effects 
on both the final-state radiation associated with the hard 
processes and the radiation induced by final-state multi- 
ple scattering in the medium. 



II. PARTON EQUILIBRATION AT RHIC 
A. Basic Equations 



densities is governed by the master equations, 

dM^) = /VR 3 (1 -A s )- 2/V R 2 (l-^), (4) 

a 

9M^) = (5) 

where i?2 = \o^n gi R3 = \ozn g , o% and 03 are ther- 
mally averaged velocity weighted cross sections, a 2 = 
(a(gg — > qq)v) and 03 = (<j{gg —> ggg)v). It has 
been calculated that R 2 « 0.064jV/a 2 A ff Trn(7.5/A 9 ) 2 
and i? 3 = 2.1a 2 T(2A s - A 2 ) 1 / 2 in Refs. HUH. 

Using Bosc and Fcrmi-Dirac momentum distributions, 
the gluon and quark densities are related to temperature 
as 



We consider here a thermal equilibrated, but chemical 
non-equilibrated system, and assume that the parton dis- 
tribution can be approximated by thermal phase space 
distributions with non-equilibrium fugacitics Ai which 
gives the measure of derivation from chemical equilib- 
rium, 



f(k;T,X i ) = X i (e^- k ±X i ) 



(1) 



where (3 is the inverse temperature and u is the four- 
velocity of the local moving reference frame, i = g,q,q 
for gluon, quark, anti-quark, respectively. 

As discussed in Ref. @ , as an approximation the mo- 
mentum distributions can be written in the factorized 
Bose or Fcrmi-Dirac form 



f(k;T,X i ) = X i (e^ k ±l) 



(2) 



which we will adopt in most of the following calculation. 

In general, chemical reactions among partons can be 
quite complicated because of the possibility of initial 
and final-state gluon radiations. However, since inter- 
ference effects due to multiple scatterings inside a dense 
medium lead to a strong suppression of soft gluon radi- 
ation, we shall only consider processes in which a sin- 
gle addition gluon is radiated, such as gg — > ggg. But 
in order to permit approach to chemical equilibrium, 
we should also consider the reverse process. Closer in- 
spection shows that gluon radiation is dominated by the 
process gg — > ggg, because radiative processes involving 
quarks have substantially smaller cross sections in pQCD, 
and quarks arc considerably less abundant than gluons in 
the initial phase of the chemical evolution of the parton 
gas. Here we are interested in understanding the basic 
mechanisms, so we restrict our consideration to the dom- 
inant reaction mechanisms for the equilibration of each 
parton flavor. These are four processes @: 



99 <-> 999, 



gg ^ qq ■ 



(3) 



Restricting to reactions, the evolution of the parton 
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P9 = ^C(3)A 9 T 3 , 



Pq 



—((3)N f \ q T 3 , 



(6) 
(7) 



where £(3) « 1.2, Nj is the number of dynamical quark 
flavors. 

If we assume that parton scatterings are sufficiently 
rapid to maintain local thermal equilibrium, and there- 
fore we can neglect effects of viscosity due to clastic 
and inelastic scatterings, we can have the hydrodynamic 
equation, 



d^eu") + Pd^ 



0. 



(8) 



where e and P are energy density and pressure of the hot 
medium. 

In order to obtain analytical solutions, we will neglect 
the transverse expansion and consider only a purely lon- 
gitudinal expansion of the parton plasma, so Eq.(8) can 
be rewritten as 



de 
rf7 



P 



0. 



(9) 



With the additional constraint of the baryon number con- 
servation, we get 



0. 



(10) 



Once the initial conditions are obtained, the evolution 
of temperature T(r), and fugacity A s (t), A g (r) can be 
determined by solving the rate Eqs.(4), (5), (9) and (10) 
together. So the input of the initial condition plays an 
important role to investigate the effects of the evolution 
for the parton system. 



B. Initial Conditions 

1. Transverse Energy 

In order to solve the rate equations discussed in the 
last section, one has to specify the initial conditions. 
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most central events of Au-Au collisions is measured after 
the time freeze out. After freezing out hard scattering sel- 
dom happens so that there is little transverse momentum 
transfer and dEx / dy\ y= $ is stable. Here only tangential 
point is the minimum of the curve and stable. If the ini- 
tial temperature is lower than 550 MeV they have two 
points of intersection, which, however, are not stable. It 
implies that the initial temperature is 550 MeV, X g o is 
0.3 for the central events of Au-Au collisions at RHIC. 
Using the same method, we obtain that the initial tem- 
perature is 420 MeV for Bjorken expansion(T 3 r = Tq to) 
for thermal and chemical equilibrium system. 



2. Particle Multiplicities 



FIG. 1: The initial fugacity A s o dependence of transverse 
energy per unit rapidity for the most central events. The 
solid, dash-dot and dash curves show the results of dEr/dy 
for initial proper time to = 0.3 fm and initial temperature 
T = 550 MeV, 650 MeV, 750 MeV at the freezing out time, 
respectively. The area between the two dot line is the data 
dEx/dyly-o = 620 ± 33 GeV for the 5% most central events 
of Au-Au collisions from RHICfU. 



The transverse energy per unit rapidity can be ex- 
pressed as [22l ] 



dE T 
dy 



tA J drid 2 p T p 2 T coshif{k;T,\i). (11) 



where £ = y — 77, y and 77 are the rapidity and pseudora- 
pidity, respectively. 

Using the momentum distribution discussed in last sec- 
tion, the transverse energy per unit rapidity can be de- 
duced as 



dE T 
dy 



[9f{\ + Ag) + 9b\ 



t-kR 2 



drjd 2 pTPr 



xcosh^e 



p T cosh£ 



(12) 



For the anti-quark fugacity, A- = . We take 
To = 0.3 fm as in Ref. @. The QGP medium freezes out 
when the temperature reduces to about 160-170 MeV. 
With different initial temperature To, we will get dif- 
ferent parton equilibration and then different curves of 
transverse energy per unit rapidity dEx / dy as the func- 
tion of A g o at the freezing out time (4.4 fm) as shown in 
Fig.ffl 

The data from RHIC shows that at ^fs = 200A GeV 
transverse energy per unit rapidity is dEx/dyly—o = 
620 ± 33 GeV[23| for the 5% most central events of Au- 
Au collisions. It is found that curve of the transverse 
energy per unit rapidity doesn't intersect the straight 
line dEx / dy\ y —o = 620 GeV if the initial temperature is 
higher than 550 MeV, only when the initial temperature 
is To = 550 MeV, X g o = 0.3, the curve is tangential to the 
straight line. dE T /dy\ y=0 = 620±33 GeU@] for the 5% 



Secondly we demonstrate that the initial condition de- 
termined above is consistent with that from the particle 
multiplicities. 

The thermodynamic functions for a many-particle sys- 
tem for an ensemble at temperature T and fugacity A can 
be derived from the grand partition, 



9fV i~ 
InZp = - orn J ap- 



6n 2 T 



1 







(p 2 - m 2 ) 1 ' 2 1 _j_ \-l e Py/p 2 +7 



1 



1 + Ae' 3 \ / P 2 +™ 2 
for fermions and 



9bV 

lnZ B = - 0rTi I dp- 



(13) 



P 



6w 2 T J (p 2 - m 2 ) 1 / 2 ^-i e /3vV+ m2 - 1 

(14) 

for bosons, where g F and g B are the degeneracy factors 
of fermions and bosons, V is the volume of the parton 
gas. 

If we neglect the mass for light quarks and gluons, 
the momentum integral Eq. (flB"]). (TTj| can be calculated 
exactly. For light quarks and gluon the grand partition 
functions can be expressed as 



TlnZ q(q) = 9 -^{\ q + X-^)T i . 
TlnZ g = ^X g T 4 . 



(15) 
(16) 



In order to get the expressions for the entropy density 
and particle number density, we recall the first law of 
thermodynamics 

E = F(V, T, fi) + TS(V, T, /x) + ^N((V, T, fi) , (17) 

where F, T, S, \x, N are the free energy, temperature, en- 
tropy, chemical potential and particle number, respec- 
tively. 

We can evaluate the free energy and the average value 
of particle number for the grand-canonical partition func- 
tion. With the thermodynamics Eq. tfTTj) . we obtain the 
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entropy density in the chemical non-equilibrated as 



h 12 , 



16. 



-lnX q (X q - \~ L )^T J - (ln\ g )\g—T . (18) 

In hydrodynamics, the relation between entropy S and 
particle number TV is S/N = 4. In a thermal equili- 
brated but chemical non-equilibrated system here, with 
considering the thermodynamics Eq. (fT7f and the grand 
canonical partition functions Eqs. (fT5|) and (|16|) . we can 
deduce 



S_ 

N 



E + PV -nN 
TN 



dNi/dy/TirR 2 



4-lnXi. (19) 



where Si and dNi/dy are the entropy density and particle 
number distribution for quarks, anti-quarks, and gluons. 
So we have 



dN/dy 



^dN^dy 

i 

_z^L_ [ g i(T3 _ (1 „ A()A , 



12 



T 3 ] 



t-kR 2 



™& ,64 3 
'4-/n(A s ) [ ^ 2 9 



k[Sa- 1 t3- ( /„a-^)a-^t3] 



16, 



(lnX g )X g — T% (20) 

7T 



where R is the transverse radius of the medium and taken 
as 6.5 fm, which is approximately the size of radius of 
Au. 

As discussed in the above subsection, with different 
initial conditions To, A g o and A g o, we obtain the differ- 
ent particle multiplicities dN/dy as shown in Fig. [21 
The data from RHIC shows that particle m ultip licities 
dN/dy = 1500 at y/s = 200A GeV from Ref. [if. Com- 
paring our calculations with the data from RHIC, we can 
conclude that the initial conditions which we get here are 
consistent with that from transverse energy. 
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FIG. 2: The initial fugacity A 9 o dependence of particle multi- 
plicities. The solid, dash-dot and dash curves show the results 
for initial proper time to = 0.3 fm and initial temperature 
To = 550 MeV, 650 MeV, 750 MeV, respectively. The dot 
line indicate dN/dy = 1500 from Ref. H|. 




FIG. 3: Time evolution of the temperature T in chemical 
non-equilibrated medium and in Bjorken expanding medium 
in Au+Au collisions for 200 GeV/nucleon at RHIC. 



C. Parton Equilibration 



With the initial conditions obtained above and rate 
Eqs. (4), (5), energy conservation Eq.(9) and baryon num- 
ber conservation Eq.(10), we can obtain the parton evo- 
lution in a chemical non-equilibrated system. The evo- 
lution of temperature and fugacity is shown in Fig. [3] 
and Fig. [4] We find that the parton gas cools faster 
than predicted in Bjorken expanding medium as shown 
in Fig. [3] because the production of additional partons 
approaching the chemical equilibrium state consumes an 
appreciable amount of energy. Fig. [4] shows that fugac- 
ity of gluons A g and quarks X q increase with increasing 
the proper time r, until the medium freezes out it hasn't 
been totally chemical equilibrated. 



III. OPACITY IN A CHEMICAL 
NON-EQUILIBRATED QGP MEDIUM 

Parton evolution will lead to the evolution of the Deby 
screening mass and parton number density. It will affect 
the cross section between the jet and the medium par- 
tons, and the opacity when the jet propagates through 
the QGP medium. 

The Debye screening mass (Id is generated by medium 
effects [Hj]. Using Bose and Fermi equilibrium distribu- 
tions in the above section, one gives the Debye screening 
mass, 
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FIG. 4: Time evolution of the fugacities X g and X q of gluons 
and quarks in the Au+Au collisions for 200 GeV/nucleon at 
RHIC. 
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FIG. 6: Time evolution of the mean free-path for gluon in 
chemical non-equilibrated medium and in Bjorken expand- 
ing medium in the Au+Au collisions for 200 GeV/nucleon at 
RHIC. 
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FIG. 5: Time evolution of the Debye screening mass in 
chemical non-equilibrated medium and in Bjorken expand- 
ing medium in the Au+Au collisions for 200 GeV/nucleon at 
RHIC. 



W r°° 

/ kf(k)dk 

Jo 

t-TTa s T 2 (\ q +2\ g ). 



(21) 
(22) 



Using the parton evolution in the above section, we can 
obtain how y, 2 D evolves with the proper time in a chemical 
non-equilibrated QGP medium and in Bjorken expanding 
medium as shown in Fig. [5j We can see that the Debye 
screening mass in the chemical non-equilibrated medium 
is less than that in the Bjorken expanding medium be- 
cause of the less fugacity in the chemical non-equilibrated 
medium. 

The cross sections for the processes, such as gg — > gg 
and gq — > gq, can be calculated from the pQCD, the 



leading order elastic scattering cross sections can be ex- 
pressed as 



9na 2 s 



27T0j 



(23) 



where /i£> is the Debye screening mass generated from the 
medium effects with considering both quark and gluon 
contributions [l9| . 

Hence, the mean free-path for a gluon l g is [25[ 



c 1 = 



a ggPg + a ggPg 
Igc(3)A s T= 



9o 



+ s C ( 3 ) NfXqT - s . (24) 



From the above results in the last section of parton 
equilibration, one can obtain the proper time depen- 
dence of the mean free-path as shown in Fig. |6l the 
inverse mean free-path decreases with increasing proper 
time, the mean free-path in the chemical non-equilibrated 
medium is larger than that in Bjorken expanding medium 
due to less Debye screening mass and larger temperature 
in the chemical non-equilibrated medium. 

We assume that jet travel the distance of L at speed 
of light in the chemical non-equilibrated QGP medium, 
the opacity can be written as 



X 



t + L 



(25) 



Using the result from parton equilibration in above 
section, the propagating distance L dependence of the 
opacity is shown in FigJT] We can see that the opac- 
ity increases with increasing the propagating distance L. 
The opacity at the initial temperature T=550 MeV in 
the static medium, which means that the medium al- 
ways stay at the initial temperature, is larger than that 
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T=550Mey- ' 




static medi«ri|i 
T=3fl"0MeV 



FIG. 7: Propagating distance dependence of the opacity in 
the Au+Au collisions for 200 GeV/nucleon at RHIC. The 
dotted curve is the opacity in a static medium for tempera- 
ture T=550 MeV, and the dotted dash curve is the opacity in 
a static medium at temperature T=300 MeV. The solid and 
dash curve indicate that in the chemical non-equilibrated ex- 
panding medium for Tb=550 MeV and the Bjorken expanding 
medium for To=420 MeV, respectively. 



in the chemical non-equilibrated and Bjorken expand- 
ing medium because of the decrease of the temperature 
and density in the expanding medium. In the chemical 
non-equilibrated expanding QGP medium, the opacity is 
larger than that in Bjorken expanding medium because of 
the higher temperature and the different ratio of quarks 
and gluons in the two mediums. One also finds that the 
opacity in the chemical non-equilibrated medium with 
the initial temperature T=550 MeV is larger than that 
in a static medium at the temperature T=300 MeV as 
considered in Ref. Hal. 



IV. PARTON ENERGY LOSS WITH DETAILED 
BALANCE 

The theoretical investigation of jet quenching are 
mostly carried out in a static medium and is found that 
the radiative energy loss is proportional to the square 
of propagating distance because of non-Abclian effect. 
However, the QGP medium is an expanding medium with 
evolution instead of being static. In this section, we will 
compare the difference of propagating distance L depen- 
dence in the different medium and the ratio of effective 
parton energy loss with and without absorption to study 
the evolution effects on jet quenching. 

In the leading-log approximation the final-state radia- 
tion amplitude of a quark can be factorized from the hard 
scattering in an axial gauge. Taking into account of both 
stimulated emission and thermal absorption in the chem- 
ical non-equilibrated expanding hot medium, one obtain 
the probability of gluon radiation with energy uj to the 



dui 



2tt 



dz 
z 



dk 2 , 
k 



, [f g (zE)S(w + zE) 



+ (1 + f g {zE))5{u - zE)6{\ - z))P gq {-), (26) 

where Cf is the Casimir of the quark jet, f g {k) is the 
gluon distribution which has been shown in Eq.(2) and 
the splitting function P gq {z) = P(z)/z = [1 + (1 — z) 2 ]/z 
for q — > gq. The first term corresponds to thermal ab- 
sorption and the second term corresponds to gluon emis- 
sion with the Bosc-Einstcin enhancement factor. 

Assuming the scale of the hard scattering as Q 2 = 
AE 2 and considering gluon radiation outside a cone with 
\kj_\ > (id> one has the kinematic limits of the gluon's 
transverse momentum [19j|, 



A<k 2 ±mm <4\ui\(E-u). 



(27) 



Subtracting the gluon radiation spectrum in the vac- 
uum, one then obtains the energy loss due to final absorp- 
tion and stimulated emission when the jet passes through 
the chemical non-equilibrated hot medium over the dis- 
tance L, 



AE 



(0) 



duiu>( 



dP (0) dP (0) 



L 

T a 

o> s CfL 

2ttL J aZ J k'^ 
+P(z)f g (zE)9(l-z)} 
Tra s C F L{T) 2 AE{T) 



=o) 



duo duj 
rlk 2 

dz I ^\-P{- z )f g ( z E) 



3L a E 



C'(2) 1 



(28) 



where 7 E « 0.5772, C'(2) ~ -0.9376, L is the thickness 
of the QGP medium, which we estimate as the diameter 
of the Au nucleon, and (T) is the average temperature 
during the time the jet propagates through the medium, 



< 



T>=lj\{r)dr. 



(29) 



As shown in Fig. [51 we see that the energy loss with- 
out rescattering increases with increasing L at fixed 
E = 5 GeV , the energy gain without rescattering in the 
chemical non-equilibrated expanding medium is larger 
than that in the Bjorken expanding medium because of 
the higher temperature in the chemical non-equilibrated 
medium. 

During the propagation of the jet after its produc- 
tion, it will suffer multiple scatte ring with targets in the 
medium. It was shown by GLV [16j that the higher or- 
der opacity corrections contribute little to the radiative 
energy loss, so we will focus on the stimulated emis- 
sion and thermal absorption associated with rescatter- 
ing with considering parton equilibration in a chemical 
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where wi = (k± - q±) 2 /2w, the factor (1 - e luliyi0 ) re- 
flects the destructive interference arising from the non- 
Abelian LPM effect. Averaging over the longitudinal tar- 
get profile is defined as <•••>= J dyp(y) ■ ■ ■ . The tar- 
get parton number distribution along the jet direction 
in the QGP medium is assumed to be an exponential 
form p(y) = 2exp(—2y/L)/L. |u(q_i_)| 2 is the normalized 
distribution of momentum transfer from the scattering 
centers, 



Kqj 
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II 2 
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d 2 (Tel 
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d 2 q± 


7T(qi- 


1 




1 


/4 




max + A*Z3 



(34) 



FIG. 8: Propagating distance dependence of the energy gain 
without rescattering in chemical non-equilibrated medium 
and in Bjorken expanding medium. 



non-equilibrated hot medium. We use the model which 
is initially proposed by Gyulassy-Wang(GW) to describe 
the interaction between the jet and target partons by a 
static color-screened Yukawa potential (2f|. Assuming 
that a parton is produced at y = {yo,yo±) inside the 
medium with yo being the longitudinal coordinate, the 
Yukawa potential is 



The propagating distance dependence of the energy 
loss and energy gain for stimulated emission and ther- 
mal absorption is shown in Fig. |9] and Fig. [10] at fixed 
jet energy E = 5 GeV. 

In the limit of EL >> 1 and E >> po, the approxi- 
mate asymptotic behavior of the energy loss of stimulated 
emission hj and energy gain of thermal absorption [l9| 
in a static medium are 



AE, 



(i) 

rad 



a s C Ff i 2 D L 2 n 2E 
'In- 



Al g E 



p 2 D L 



0.048] , (36) 



V„ 
v(q n ) 



2ir5(q Q )v(q n )e-^-yH an (j)t an (n), (30) 

(31) 



9 i 



where q n is the momentum transfer from a target parton 
n at the position y n = (y n ,y n ±), t a „(j) and t a „(n) are 
the color matrices for the jet and target parton. 

Here we consider the contributions to the first order 
in the opacity expansion. By including the interfer- 
ence between the process of the rescattering and non- 
rescattering, we obtain the energy loss for stimulated 
emission AE^J d and energy gain for thermal absorption 
Ae[\\ to the first order as 



abs 



AE, 



(i) 



a s C F E r a+L 



v(q±) 



dr 

qj 



dk\ 



dz 

P{z){a gg p g 



AE 



(i) 

abs 



(k± - q±_ 
< Re(l - e iwiyi °) >0(l-z), 
a s C F E [ Ta+L 



d 2 q A 



VgqPq) 



(32) 



dr 

qj 



dz 



dk\ 



d 2 q, 



7fg{zE){aggp g + a gq p q ) 



I I 2 n 

(k± - qj_) z 
x [-P(-z)) < Re(l - e lUlVl °) > 
+P(z) < Re(l - e iuJiyi °) > 9(1 - z)]. (33) 



AE 



(i) 



E 



ira s C F LT 2 p 2 D L 

;[ln-^ 1- 



T 



-lE- 



6C (2)! 



(37) 



As shown in Eq. (|36p . the analytic approximation of the 
energy loss by stimulated emission in a static medium 
is proportional to L 2 ln(j-). We fit the curve from 
our numerical calculation in a chemical non-equilibrated 
medium, it is shown as in Fig. [9] that the energy loss 
by stimulated emission is proportional to L approxi- 
mately by taking into account the parton evolution of 
the chemical non-equilibrated medium and Bjorken ex- 
panding medium. In a static medium it is shown that the 
energy gain AE^ S in Eq. (|37[) , from thermal absorption 
is linear distance dependence if assuming p 2 D L/T >> 1, 
and become a quadratic distance dependence if assuming 
p 2 D L/T « 1. In our case for a chemical non-equilibrated 
medium here, from the curve fitting it is shown as in Fig. 
[TU] that the energy gain keep the linear distance depen- 
dence as p 2 D L/T >> 1, but the distance dependence be- 
comes L 2 as p 2 D L/T « 1 instead of L 2 -dependcncc 
in static medium case and the proportionality factor is 
much smaller than that in the static medium. 

As shown in Fig. [H] the energy loss without ther- 
mal absorption in a chemical non-equilibrated medium 
with initial temperature Tq = 550 MeV, initial fugacity 
Xgo = 0.3 is larger than that in a static medium with 
T = 300 MeV, pi D =0.5 GeV and the mean free path 
of the jet l g = 1 fm. But it is less than that in the 
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FIG. 9: Propagating distance dependence of the energy loss 
to the first order opacity. For static medium case, T — 300 
MeV, fiD = 0.5 GeV and l g = 1 fm. For chemical non- 
equilibrate medium case, To = 550 MeV , X g o — 0.3 and to — 
0.3 fm. For the Bjorken expanding case with two different 
initial temperature, T = 420 MeV and T = 550 MeV, 
to = 0.3 fm. 
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FIG. 10: Propagating distance dependence of the energy gain 
to the first order. For static medium case, T = 300 MeV , 
)Xd = 0.5 GeV and l g = I fm. For chemical non-equilibrate 
medium case, Tb = 550 MeV, X a o = 0.3 and to = 0.3 fm. 
For the Bjorken expanding case, To = 420 MeV, to = 0.3 
fm. 



Bjorken expanding medium with T = 420 MeV . Al- 
though the opacity in the Bjorken expanding medium is 
less than that in the chemical non-equilibrated expanding 
medium, the Debye screening mass is larger, this leads to 
larger q_ max and larger energy loss per collision. We also 
calculated the energy loss without thermal absorption in 
the Bjorken medium with the same initial temperature 



FIG. 11: The ratio of effective parton energy loss with and 
without absorption as a function of parton energy E and prop- 
agating distance L. 
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FIG. 12: The ratio of effective parton energy loss with and 
without absorption as a function of propagating distance L 
at fixed E — 3 GeV in a static medium and chemical non- 
equilibrated hot medium cases. 



To = 550 MeV as in chemical non-equilibrated medium 
to study difference between Bjorken expansion and chem- 
ical non-equilibrated expansion medium. It is shown 
that the energy loss without thermal absorption in the 
Bjorken medium is larger than that in the chemical non- 
equilibrated medium due to that temperature decreases 
more rapidly in the chemical non-equilibrated medium 
than the Bjorken expanding medium. As shown in Fig. 
1101 the energy gain in the chemical non-equilibrated 
medium is a bit larger than that in the Bjorken expanding 
medium with To = 420 MeV due to the larger tempera- 
ture in the chemical non-equilibrated medium. 
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Fig. [TT] shows the ratio between parton energy loss 
with thermal absorption AE = AEi ra( i + A-Eia&a + 
AEo a b s and the radiative energy loss AEi ra( i with- 
out thermal absorption as functions of parton energy 
value E and propagating distance L in a chemical non- 
equilibrated medium. It is shown that the energy gain 
via absorption is important for intermediate parton en- 
ergy and small propagating distance. The absorption 
can be neglected either for asymptotically large parton 
energy or large values of L. The comparison between 
the ratios in a static medium and in a chemical non- 
equilibrated hot medium at fixed parton energy E = 3 
GeV is shown in FigfT2l We can see that the ratio of the 
calculated parton energy loss with and without thermal 
absorption in a static medium(T = 550 MeV, fj,D = 1 
GeV) considered in Ref. is less than that in a chem- 
ical non-equilibrated hot medium. This result implies 
that the energy gain via absorption is more important in 
a static medium than that in a chemical non-equilibrated 
medium at large distance L and intermediate jet energy 
region. 

V. CONCLUSION 

In summary, we have determined the initial conditions 
for central Au-Au collisions at RHIC. It is shown that 
the initial temperature is 550MeV, the initial fugacity 
X g for gluons is 0.3 and \ q for quarks is 0.06 for the 
chemical non-equilibrated expanding medium. The ini- 
tial temperature is 420 M eV for the Bjorken expanding 
medium. We obtain that partons cool faster than pre- 
dicted by Bjorken with the proper time. This chemi- 
cal equilibration in a hot medium lead to that the De- 
bye mass and the inverse mean free-path decrease with 
increasing the proper time, and the opacity for a en- 
ergetic parton in a chemical non-equilibrated medium 
increases with incraesing the propagating distance L. 
The Debye mass and the mean free-path in the Bjorken 
expanding medium is larger than that in the chemical 



non-equilibrated expanding medium. We found that the 
opacity in the chemical non-equilibrated medium is larger 
than that in the static medium at temperature 300 MeV 
and that in the Bjorken expanding medium, but is less 
than that in a static medium at T = 550 MeV . 

The parton evolution affect the energy loss with de- 
tailed balance. The propagating distance L depen- 
dence of parton energy loss with detailed balance is 
determined. The net effect of stimulated gluon emis- 
sion and thermal absorption are considered for an en- 
ergetic parton propagating inside the QGP medium. It 
is shown that, by taking into account the evolution of 
the temperature and fugacity, the energy loss from stim- 
ulated emission is proportional to L in the chemical non- 
equilibrated medium and Bjorken expanding medium 
rather than L 2 -dcpendence on the propagating distance 
in the static medium. The energy loss in the chemi- 
cal non-equilibrated medium is a bit less than that in 
Bjorken expanding medium. Because of the higher tem- 
perature in the chemical non-equilibrated medium, the 
energy gain without rescattering to the first order opac- 
ity in the chemical non-equilibrated medium are larger 
than those in the Bjorken medium. The relative reduc- 
tion of the parton energy loss is significant at intermedi- 
ate parton energy E and small propagating distance L in 
the chemical non-equilibrium medium. The evolution of 
the medium modifies the jet energy loss in the interme- 
diate energy region and affect the shape of suppression 
intermediate high Pt hadrons spectrum. 
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